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Abstract. In this announcement we consider an eigenvalue problem which 
arises in the study of rectangular membranes. The mathematical model is 
an elliptic equation, in potential form, with Dirichlet boundary conditions. 
We have shown that the potential is uniquely determined, up to an additive 
constant, by a subset of the nodal lines of the eigenfunctions. A formula is 
given which, when the additive constant is fixed, yields an approximation to 
the potential at a dense set of points. An estimate is presented for the error 
made by the formula. 

Introduction 

In this summary we consider the Dirichlet eigenvalue problem for the operator 
(— A +q) on a rectangle. This problem arises in the study of rectangular, vibrating 
membranes. The goal here is to solve the inverse problem: find q from the nodal 
line positions of the eigenfunctions. We show that, for almost all rectangles and for 
sufficiently smooth q , a potential q whose integral average is zero can be uniquely 
determined from a subset of the nodal lines of the eigenfunctions. 

The theorems, which we present in this paper, extend the one-dimensional re- 
sults of McLaughlin and Hald [8], [4], [5], where inverse nodal problems are de- 
fined. There the authors showed that the one-dimensional potential for the Sturm- 
Liouville problem with Dirichlet or mixed boundary conditions is uniquely deter- 
mined, up to an additive constant, by a dense set of nodes. The analysis of the 
Sturm-Liouville problem depends on perturbation theory. The proofs of the per- 
turbation results are reasonably straightforward. This is possible because the dif- 
ference between consecutive eigenvalues increases as the order of the eigenvalues 
increases. The perturbation results provide asymptotic forms for large eigenvalues, 
the corresponding eigenfunctions and nodal positions. The uniqueness theorem 
follows. 

For the two-dimensional problem we also require asymptotic forms for the eigen- 
values and eigenfunctions and approximate location of nodal lines. The difference 
between the one-dimensional and the higher dimensional cases is that in the higher 
dimensional case the eigenvalues are not well separated. Even for those rectangles 
for which all eigenvalues are distinct, arbitrarily large eigenvalues can be arbitrarily 
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close together. These small differences can become small divisors in an asymptotic 
expansion; see Kato [7]. To handle this, we start with the q — case and show 
there is a dense, well-defined set of rectangles for which the eigenvalues are well 
separated. The rectangles are selected so that the square of the ratio of the sides 
is not well approximated by rational numbers. Our condition is almost the same 
as that given by Moser [9]. For each acceptable rectangle, we choose "good" eigen- 
value, eigenfunction pairs that have the following properties. These eigenvalues are 
at least a given distance from their nearest neighbors and a large distance from a 
selected set of neighbors. The corresponding eigenfunctions have nearly the same 
number of oscillations in the x and y directions. Using number theoretic and 
geometric results, it can be shown that almost all eigenvalue, eigenfunction pairs 
satisfy these properties. 

Using the q = case as a base problem, the application of perturbation the- 
ory yields leading order terms for "good" eigenvalue, eigenfunction pairs for the 
q =/= case. Explicit estimates are obtained for remainder terms. Our technique is 
strongly influenced by perturbation results for the two- and three-dimensional peri- 
odic problem by Feldman, Knorrer, and Trubowitz [3] and Friedlander [2] . Instead 
of L 2 estimates, we need and derive L°° estimates for the difference between the 
eigenfunctions and their leading order terms; we require the L°° bound so that we 
can determine the approximate location of the nodal lines for the q ^ case. 

Having established asymptotic forms for almost all the eigenfunctions and having 
established the L°° bound for the difference between the eigenfunctions and their 
leading order terms, we encounter the second major difficulty. A nodal domain for 
an eigenfunction for the q ^ case is not necessarily a small perturbation of a 
nodal domain for the corresponding eigenfunction for the q — case. This follows 
from the fact that the nodal line intersections for the q = case are not usually 
present in the nodal line patterns for the q ^ case. The fact that nodal line 
intersections are not stable under perturbation has been shown by Uhlcnbeck; see 
[11], [12]. Figure 1 on page 245 shows typical nodal domains for the corresponding 
eigenfunctions for the q = case and the q ^ case. 

In our analysis we select approximate nodal domains which are subsets of nodal 
domains for the q ^ problem. The approximate nodal domains are small per- 
turbations of nodal domains for the q = problem. The eigenfunction for the 
q 7^ case is either zero or nearly zero on the boundary of the approximate nodal 
domain. 

In all of the above, we address properties of the forward problem. Only then can 
we address the inverse problem. 

We solve the inverse nodal problem. A formula is given for approximating q at 
a dense set of points using only a subset of the nodal line positions. We show that, 
at each point, q can be approximated by the difference of two eigenvalues. One 
eigenvalue is an eigenvalue for the q = Dirichlet problem defined on the rectangle. 
The other eigenvalue is the first eigenvalue for the q = Dirichlet problem defined 
on an approximate nodal domain. An estimate is given for the error in making the 
approximation. A uniqueness theorem follows. The proofs, which are presented in 
[6] , use the Raylcigh-Ritz variational formulation of the eigenvalue problem and the 
perturbation results described above. 

In the body of this paper we present statements of our key results. Here we do 
not give the sharpest possible results but give a presentation that makes the ideas 
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more accessible to the reader. The complete set of results together with the proofs 
are contained in [6]. 



Main results 

We begin with the problem 

(1) — Au + qu = Xu , (x , y) e R, 

u = 0, (x , y) e dR , 

where R = [0 , it /a] x [0 , it] , a > 1 and q e Cfi° (R) . Our goal is to solve the 
inverse problem: Find q from the positions of the nodal lines. To attain this goal, 
we must first study the forward problem. When q = , that is, when we have the 
problem, 

— Au = Xu, (x , y) e R, 
u = , (x , y) € dR , 

we can label the eigenvalue, eigenfunction pairs as 

Aq,o = \a\ 2 = \(an , m)| 2 = a 2 n 2 + m 2 , 

u a o = sin anx sin my , 

respectively. Here a is in the two-dimensional lattice 

L(a) = {a = (an , m)| n,m=l,2,3,...}. 

We must first establish a perturbation result for the perturbation from q = 
to ij ^ 0. We can do this for almost all of the eigenvalues {X a a} a eL when the 
parameter a , which defines the rectangle R , is well chosen. Specifically we choose 
a 2 so that it is not well approximated by rational numbers; i.e. a is in the set 

V = I a > 1 1 there exists < S < — and k > , 

such that for all p, q > : la 2 — -I > ^ x > . 

q q z+a I 

It can be shown that in any subinterval, J = (1 , a ) , mcas (J\V) = 0. Further 
we can write down specific irrational numbers which are in V . For example we 
can show that \J e — 1 e ; and it follows from a theorem of Roth [10] that if a 2 
is an irrational, algebraic number, then a . 
Our perturbation result is 

Theorem 1. Let a belong to V . Let q e Cfi°(R) have mean value 0. There 
exists an exceptional set of lattice points M(a) C L(a) depending on a and an 
upper bound on the derivatives of q of order no more than 100 such that 
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a) M(a) has density in L(a) in the sense that 



Um #{a e L(a)\M(a) \ \a\ < r} = ^ 
r->oo #{a e L(a) | |a| < r} 

b) For every a G L(a)\M(a) , there is a unique eigenvalue of the variable coef- 
ficient problem (1) which we denote as \ aq satisfying 

\x aq -x a o \< i«r 15/s . 

A suitable multiple of the corresponding eigenfunction satisfies 



E(qu a0 , u [i0 ) 
U2_|fl|2 U P° 



< Va\c 



-15/8 



Remark. Our techniques for the proof of Theorem 1 are strongly influenced by [1] , 
[2], [3], [7], [13]. 

Remark. We note that there are two properties of the L°° bound for the eigcn- 
functions that are important to us. One is that we choose "good" a so that we 
have 



E(qu a0 , upo) 



a • 



l/?l 5 



-15/16 



We use this to establish the position of the nodal lines for u aq . The second property 
is that near the intersections of the nodal lines of u a0 the bound 



< y/a\a\ 



-15/8 



holds; this is used to establish the error bounds for the formula for the potential 
q. 

We have established the L°° bound for u aq — u a o (see the Remark above) so 
that we can determine the approximate location of nodal lines and hence define 
approximate nodal domains for u aq . We first exhibit a nodal domain for u a o , 
a = (an , to) ; it is 



On 



(x, y) 



(m + !)tt 



an 



< 



2an 



(mi + i)7r 



TO 



< 



TT 

2n~i 



for some n\ =0,1 



1 , TOi =0,1, 



1 , and note that a nodal 



domain of u aq with aeL\M is a connected component of {(x, y)\u aq j^0} . Now 
we have a difficulty because a nodal domain of u aq may be a large perturbation of 
a nodal domain for u a o ■ Figure 1 shows the nodal domains for u Q n and u aq for 
a typical case. 

Note that the nodal line intersections for u q0 are not present in the nodal line 
pattern for [11], [12] for the proof that nodal line intersections are unstable 

under perturbation. 

In Figure 1, the "natural" perturbation of fi is the much larger nodal domain 
£1 . We will select a subdomain of f2 as our approximate nodal domain. To do 
this more precisely, let d = l/[(3.9)|a| 15 / 16 ] , d\ — d/(an) , d 2 = d/m, and let 
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nodal domains for u a o 

fii C Ct C Q, 2 be 
Oi = \ (x, y) 



nodal domains for u a 



(m + ±)tt 



2 



x — 



an 
(ni + !)tt 



an 



Figure 1 



2an 1 



< 2an ^ 



(mi + |)7r 



m 

(mi + ±)tt 



m 



. This allows us to choose a natural nodal 



Figure 2 exhibits the three domains Qq , 
We can then show that u aq ^= in £l\ 
domain Q associated with fig . It is the nodal domain of u aq which contains Q,\ . 
Our approximate nodal domain is a subset of O and is as follows. 

Definition. Let a e L(a)\M(a) and let Slo be a nodal domain for u a0 . Define 
the corresponding Sli C SJo C as above. Let Q be the nodal domain of u aq 
which contains f2i . Then the perturbation £1' of Qq is the connected component 
of £1 n VI2 which contains Oi . 

Figure 3 depicts a choice for O' . 

These are the main theorems of this paper. 

Theorem 2. Let q e Cq°(R) with mean value zero, and let a e L\M . Let Oo 
be a nodal domain of u a0 , and let Q' be the approximate nodal domain of u aq 



Figure 2 
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Figure 3. Construction of 0' . 
defined above. Let A 10 (^'J be the first eigenvalue of 

-Au = Xu , (x , y) £ Cl' , 
u = 0, (x , y) £ dQ,'. 

Then, there exists (x , y) <G with 

\q(x,y)-[\a\ 2 -\ 1 , (tt)}\ < l«r ?/4 - 

Because M has density in L , the set of points (x , y) for all nodal domains of 
u a o , a G L\M , is dense in R . This proves 

Theorem 3. Let q e Cfi°(R) with mean value zero and let a G L\M . Then q is 
uniquely determined by the nodal lines of the eigenf unctions for 

(- A + q)u = Xu , (x , y) e R, 

u = , (x , y) <G dR. 

Remark. We can improve on Theorem 3 showing that q is uniquely determined by 
a subset of the nodal lines. To do this we use the property that, when a e L\M , 
an and m are comparable. From this it follows that q is uniquely determined by 
the nodal lines of u a k q for any set a k £ L\M , k = 1 , 2 , ... where | a k |— > oo as 
k — > oo . 
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